Weak solutions to lubrication systems describing 
the evolution of bilayer thin films 

Sebastian Jachalski* Georgy Kitavtsev^ Roman Taranets* 

November 12, 2012 



The existence of global nonnegative weak solutions is proved for cou- 
pled one-dimensional lubrication systems that describe the evolution of 
nanoscopic bilayer thin polymer films that take account of Navier-slip or 
no-slip conditions at both liquid-liquid and liquid-solid interfaces. In ad- 
dition, in the presence of attractive van der Waals and repulsive Born in- 
termolecular interactions existence of positive smooth solutions is shown. 

1 Introduction 

During the last decades lubrication theory was successfully applied to modeling 
of dewetting processes in micro and nanoscopic liquid films on a solid polymer 
substrates see e.g. [H [TSl [T7] to name a few, for a review we refer to [5] and 
references therein. A typical closed-form one-dimensional lubrication equation 
derived from the underlying equations for conservation of mass and momentum, 
together with boundary conditions for the tangential and normal stresses, as 
well as the kinematic condition at the free boundary, impermeability and a slip 
condition at the liquid-solid interface has the form: 



where function h(x, t) denotes the height profile for the free surface of the film. 
The mobility function has the form M(h) = h 3 or M{h) = bh 2 for the no-slip 
or Navier-slip conditions considered at the solid-liquid interface, respectively, 
where b > denotes the slip-length parameter. 

Recently, this model was generalized to a coupled lubrication system de- 
scribing evolution of a layered system of two viscous, immiscible, nanoscopic 
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Abstract 




(1.1) 



Newtonian fluids evolving on a solid substrate [TJ Q3] and subsequently anal- 
ysed in [2] [13j [TU [16j [18] . The latter system can be stated in the form: 

u t = -d x {MudxPi + M 12 d x p 2 ) , ^ 
v t = -d x {M 12 d x px + M 22 d x p 2 ) , 

where u(x, t) and v(x, t) denote the height of the lower liquid and the difference 
between the heights of the upper and lower liquid, respectively (see Fig. 1). 
The pressures pi(x, t) and p 2 (x, t) are defined as 

p 2 = dlu + d 2 x v-U 2 (v), 

where d x u and d x v are linearised surface tension terms and potentials Hi(u) 
and H 2 (v) describe the intermolecular interactions of the bottom liquid with 
the solid surface and of two liquids with each other, respectively. The influence 
of intermolecular interactions is typically due to the competition between long- 
range attractive van der Waals and short-range Born repulsive intermolecular 
forces, see [10l [17]. In this article we consider two case: the absence of inter- 
molecular interactions, i.e. Ilfc(s) = for k — 1, 2 and the case when both van 
der Waals and Born intermolecular forces are presented in the form 

n fc( s )^7T-^' (™< m < 7i,72 «1) (1.4) 

A typical choice for (n, m) is (3, 12) corresponding to the standard Lennard- 
Joncs potential. 

As in the case of the single layer lubrication equation (jl.lj) the form of 
mobility matrix 



M (u, v) 



f M u (u,v) M 12 (u,v)^ 
M 21 (u,v) M 22 {u,v), 



depends on the slip conditions considered at the liquid-liquid and liquid-solid 
interfaces. In the case of the no-slip at the both interfaces, it has the form 




1 2 
2 U V 



M = - . (1.5) 



3 u 



The model parameters a = a±/a 2 and p = pi/fJ-2 in (| 1 -3[) and (jl.5l) are positive 
constants which denote the ratios of surface tensions and viscosities, respectively. 
Recently, the lubrication system for the case of Navier-slip conditions considered 
at both liquid- liquid and liquid-solid interfaces was derived in [HI IT8] . The 
corresponding mobility matrix can be stated in the form 
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Figure 1: Sketch 



where b\ > and b > denote the slip lengths at the solid-liquid and liquid- 
liquid interfaces, respectively. We will discuss the origin of the different mobili- 
ties in section 5. Note that rescaling time by b\ and introducing the parameter 
a := jpOu + 1) > the latter matrix can be written in the form 

, / u 2 uv \ 

M=-\ . (1.6) 

f* \uv (l + a)v 2 J 

In this study we consider the system (11.21) on a time-space domain Qt = 
ft x (0, T) where O = (0, 1) with boundary conditions 

u x = u xxx — v x = v xxx =0 for x e dQ (1.7) 

and the initial functions 

u(x,0) = u (x) > 0, v(x,0) = v Q (x) > u^vell 1 ^). (1.8) 

We consider the system (|1.2p both with the mobility matrixes (|1.5I) and p.6[) . 
The system (|1.2p can be also generalised to incorporate presence surfactants or 
temperature-gradient-caused Marangoni flows (see e.g. [71 1^1 HB]). 

Starting from the seminal work of Bernis and Friedman [5] existence the- 
ory of nonnegative weak solutions for single nonlinear parabolic equations in 
the form was successfully developed (see e.g. [6] and references there 

in). Note that the system (|1.2j) inherits from the one-layer lubrication equation 
(II. 1|) the high order and degeneracy as one of the fluid heights u or v goes to 
zero. In contrast to there are only few analytical results known about 

the system (|1.2[) . The structure of its stationary solutions were considered in 
[131 HE]- Existence of nonnegative weak solutions to (|1.2p in the no-slip case, 
i.e. with mobility matrix given by (ll.5[) was shown recently in [2] using a finite- 
element approximation under a strong assumption on presence of intermolecular 
potentials of the form (|1.4p between liquid films and between each film and the 
substrate as well. 

In this article, we show existence of nonnegative weak solutions to (|1.2p with 
(ll.7p - (ll.8p for both no-slip f|l .5[) and Navier-slip cases (|1.6p in the absence of 
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intermolecular forces. In turn in the presence of intermolecular forces between 
liquid films and just between bottom liquid and the substrate as in (|1.3[) we 
show that the observed weak solution becomes positive and smooth. 

In our approach we extend ideas introduced in [3] for the single lubrication 
equation of the form (II. ip to the system (|1.2j) . However, the extension is not 
straight forward. There are new challenges since the mobility matrix degenerates 
in more than one way. Beside the case that u and v vanish simultaneously, 
also the mobility matrix M may degenerates if either u or v becomes zero 
while the other does not. All these cases have to be treated very carefully 
In section 2 we introduce the corresponding regularized version of the system 
(jl.2p for the no-slip case, with mobility matrix (jl.5p . By deriving the energy 
dissipation and corresponding a priori estimates, using the theory of uniformly 
parabolic systems (see [11]), we show global existence of smooth solutions to 
the regularised problem. Furthermore, we show that the latter converge to a 
suitably defined weak solution of the original system (|1.2|) . Notice that these 
results are independent of the presence of intermolecular forces in the equations. 

In section 3 we prove nonnegativity of thus obtained weak solutions in the 
case when intermolecular forces are absent. In the presence of intermolecu- 
lar forces as in Q1.2p the weak solutions turn out to be positive and smooth. 
Our approach for proving the nonnegativity is based on a definition of suitable 
analogs of Bernis and Friedman entropies for functions u and v and showing 
their combined dissipation. 

Moreover, in section 4 we show global existence of nonnegative weak solu- 
tions to (|1.2p in the Navier-slip case (|1.6p . 

2 Existence of weak solutions in the no-slip case 

In this section we consider the system (|1.2p without intermolecular interactions, 
i.e. with pi,p2 in (|l-3p given by 

pi = (a + l)u xx + v xx , P2 = u xx + v xx ; (2.1) 

and with the no-slip mobility matrix 

IN 3 §M>I \ 
§M>| fH 8 + |«IHv ' 

Notice that we replaced every u and v in the mobility matrix (jl.5p by their 
absolute values to ensure that the latter is positive semidefinite. We will prove 
existence of global weak solutions to (|1.2p , (|2.ip considered with boundary and 
initial conditions (|1.7|) - (|1.8|) . We begin our analysis with introduction of a 
regularised version of (|1.2p . (|1.8I) and derivation of a priori estimates for its 
solutions. 
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2.1 Regularised system and a priori estimates 

Since (jl.2p is degenerate at u — and v = 0, we approximate it by a family of 
non-degenerate equations 

Ut + ((Mn + £)Pl,z + Mi 2 P2,x)a: = . , . 

v t + (M 2lPhx + (M 22 + e) P2 , x ) x = 

where £ > is arbitrary. Note that the regularised mobility matrix is positive 
definite for all u and v. Correspondingly the system (|2.2[) is uniformly parabolic 
in Petrovskii sense (see jlll for definition). Furthermore, we approximate uq 
and Vo in the H 1 (O)-norm by C 4+a functions uq £ and vq £ satisfying (|1.7I) , 

Moe(aO > uo(x) and Vo e (x) > wo(a^) for x 6 f2, (2-3) 

and replace (|I.8[) by 

u(a:,0) = uoe(a:), u(a;,0) = v 0e (a;). (2.4) 

By [HI Theorem. 6.3, p.302] the system considered with ([T7r|) . (gH) has 
a unique local solution (u £ ,v £ ) in Q T for some small r = r(e) > 0. 

Everywhere below in this article we denote by C positive constants indepen- 
dent of e which may vary from line to line. Let us also introduce notations 

M e = M(u e ,v e ) 

Pi,s = (c + l)u e , X x + v e>xx , and p 2:£ = u e>xx + v £ , xx . (2.5) 
A priori estimates 

Let us define an energy (Lyapunov) functional for the system (|1.2[) coupled with 
PTTf]) as 

#(u E , v e ) = y [<Ju^ x + (u £ , x + v £iX ) 2 ] dx (2.6) 
h 

Indeed, differentiating the latter in time along solutions of (|I.2p with (|1.7p one 
obtains the corresponding energy equality 



1 d . 

— —E(u £ . v £ 

2dt v 



J {Mn £ p\ ex + 2M l2 ePle,xP2e,x + M 22£ p\ ex )dx 

h 

+ £ J (Ple,x + pL,x) dx = °- 



(2.7) 



Note that the second term in (|2.7[) is nonnegative since M e is positive semidef- 
inite. By the approximation properties of uq £i vq £ one has 



u 0e , x dx < (1 + T](e)) I u x dx, 

(2.8) 



v e,xdx < (! + ? 7( £ )) / v o,x d x 
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where 77(e) — > if e — > and therefore E(uo e ,VQ e ) < C holds. This together 
with (I2.7P imply the following a priori estimates: 



sup / u 2 £ x dx < C, sup / u 2 . x dx < C, 

;E(0,t)J ' *G(0,r)J 



te(o,r)J te(o,- 



(2.9) 



and 

/ f (MiuPie,x + 2M 12ePl ^ xP2 ^ x + M 22e p 2 2e>x )dx < C, (2.10) 

Qt 

e JJ (p 2 le , x +pL x )dx<C. (2.11) 

Integrating (|2.2[) in time we deduce the conservation of mass law 

u £ , x (x,t)dx = u e tX dx, / v s>x (x,t)dx = v 0e , x dx (2.12) 



for all t £ (0,t). Using this, ()2.8|) . Poincare's inequality and the Sobolev em- 
bedding theorem 

one obtains 

IK(.,0H c o 4(fi) <^ lk(.,0ll c o., (S) <C. (2.13) 

Next, we obtain uniform Holder estimates for u e and v e in time. Let us 
introduce functions 

Jl )£ = MliePie.aj + M 12e P2e,x and J 2)£ = M21 £ pi £)X + M 22e p 2s>x . 

Observe that for every t G (0, r) 

Jj 2 e dxdt 
Qt 

<C J J M U e (M lle p\ ex + 2Mi2ePle,xP2e,x + M 22eV \ ex ) dxdt 

Qt 

< C, (2.14) 
where we use Ml 2e < M lle M 22e , (j2~T0t and (l2~T3|) . Analogously, 

Ji e dxdt<C (2.15) 
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holds. Now, using (|2.14j) ~ (l2.15l) and the relations 



u e 4> t =11 J\, e 4>x and / / v E 4> t = - 

Qt Q 



>2,e<Px 



considered with the special test function <j> taken exactly as in the analogous 
proof for the single layer lubrication equation in j3J Lemma 2.1] one obtains 
that for all x € Cl and t%, t 2 in (0, r) the following holds 

\u e (x,t 2 )-u e (x,t 1 )\<C\t 2 -ti\i, 

i I 2 - 16 ) 
\v e (x,t 2 ) - v £ (x,h)\ < C\t 2 -h\s 

Conclusion 

The relations (|2.13[) and (|2.16p imply upper bounds on the C£' t a -norms of u e 
and v e in Q T , which are independent of r, e. These a priori bounds allows us 
to conclude that (u e ,v e ) can be extended step-by-step to a solution of (|2.2[) 
considered with (|1.7[> . (|2 .4[) in Qt for any positive T > (see [TTJ Theorem. 
9.3, p. 316]), and that 

the sequences {u £ } and {v e } are a uniformly bounded 
and equi-continuous families in Qt- 

2.2 Convergence to global weak solutions 

Here we show that solutions u £ , v e of the regularised system (|2.2[) converge to 
suitably defined global weak solutions of the initial system (|1.2p . By (I2.17p . 
every sequence e — > has a subsequence (for short both not labeled) such that 

u e — > u, v e — > v uniformly in Qt- (2-18) 

Note that, due to uniform bounds (I2.13[) and (12 . 16[) . any such limits u and 
v can be defined globally in time using a standard Cantor diagonal argument 
(choosing a sequence T n — > oo). 

Theorem 2.1. Any pair of functions (u,v) obtained as in (2.18)) satisfies for 
any T > the following properties: 

u,ve C^ 2 ' 1/8 ((5t), u,v e C 4 x j(P), (2.19) 
M:\P .r + M 12 p2,x, M 21 p hx + M 22 p 2 , x e L 2 (P), (2.20) 
\u\W x E L 2 (R), \v\ 3 p 2 , x e L 2 (S); (2.21) 

where_P = Q T \{{u = 0} U {v = 0} U {t = 0}), R = Q T n{v = 0}n{|u| > 0} and 
S = Qt H {u = 0} n {\v\ > 0}. Furthermore, there exists a function w € L 2 (R), 
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such that (u, v) satisfies U.ty) in the following sense: 

ucf>t+ J J (M n p ltX + M 12 p 2 , x ) 4>x (2.22) 

Qt P 

R 

v<p t + J J {M 2lPl . x + M 22 p 2 , x ) <f> x + J J ^\v\ 3 p 2 , x 4> x = (2.23) 

Qt P S 

for all <p G Lip(Qx), 4> — near t — and t = T: 

u(x, 0) = uq(x), v(x, 0) = vo(x), xe£l, (2.24) 

\\u(-,t)\\LHn) = IKIUi(n)) ||«(-,*)||ii(n) = ||wolUi(Q), (2-25) 

u x (-,t) — > u 0x and v x (-,t) — > v x strongly in L 2 (fi) as t — > 0, (2.26) 



and 



u and v satisfy \1. 7[ ) at aH points of the lateral 
boundary, where u ^ and ti ^0. 



(2.27) 



Finally, the following energy inequality holds 

E(u(;T),v(-,T))+ ff {M 2lP l x + M 22 p\ x ) <E(u ,v ). (2.28) 



p 



Proof. By the properties of the constructed solutions u £ , v e to the regularised 
system and their uniform convergence to u, v the first assertion in (|2.19l) and 
also (|2.24l) - (|2.25p follow immediately. Using (I2.10|) . we observe 



(Mi l£ p lex + M 2 2 £ p 2ex )dxdt < C - 2 J J M X 2eP\e,xP2e,xdxdt 

Qt 

^ J J \u e \ 2 \v e \pu >x p2 S , x dxdt, 



Qt 

= c 



Qt 



From this applying Young's inequality 

K\ 2 \v £ \ PlE , xP 2s,x < -^K\ 3 p 2 le . x + -KIN^L,* 

one obtains 

{M lu p\^ x + M 2 2eP2e, x )dxdt < C (2.29) 
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and therefore estimates 

\u E \*p\ E>x dxdt <C, J J \v e \ 3 pl EtX dxdt <C, J J \u e \\v E \ 2 p 2 EtX dxdt < C. 

Qt Qt 

(2.30) 

For any <f> as in (I2.22[) one has 



u e <p t dxdt + J J ((Mii E +e)px e> x + M 12s p2 e , x ) <t>xdxdt = 0, (2.31) 

Qt 

v e 4> t dxdt+ II {M 21ePle . x + (M 22e + e)p2e.x) (f> x dxdt = 0. (2.32) 



Qt Qt 

By (|2.14l) and (|2.15l) there exist Ji,J 2 € L 2 (Qt) and a subsequence as e —> 
such that 

(MuePle,x + M 12e p 2e . x ) -» Ji 

and (M 2l£ p letX + M 22s p 2e ,x) J 2 weakly in L 2 {Q T ). (2.33) 
Additionally by (l2~TTj) 

efjp^xdxdt - 0, eJJ p 2e ^ x dxdt - as £ ^ 0. 

Qt Qt 

By regularity theory of uniformly parabolic systems and the uniform Holder 
continuity of the u E and v E we deduce that C x 'l(P) and 

Ji = MhPi, x + Mi2j?2,x, J? = M 2 iPi,x + M 22 p 2tX in P. (2.34) 

Next, for a fixed <5 > define a set Ii.s = {\v\ < S < \u\}. From the 
estimates (12.30)) it follows that there exists w € L 2 (Ii t s) such that Pi e , x — ^ 
and v e p 2e . x — 1 iy weakly in L 2 {I\^) as e — >• 0. Therefore, one obtains 

^ (MnePiE.x + Mi2e_p2 £ ,^) 4> x dxdt -> ^ (^M^M + 2~l u ' 2u; ) <t>xdxdt 

Il,S Il,S 

(2.35) 

as e — ?> 0. On the other hand, one has the estimate 
[[ (M lle pi E<x + M 12E p 2eiX )4> x dxdt 

\u\<5 

< C J J M XXe (M Ue p 2 EtX + 2M 12ePi 

e,xP2e,x + M 22e p\ E x ) dxdt 

\\u\<6 

< CS 3/2 (2.36) 
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Let us decompose the second term in (|2.31[) as follows 



(Mu e pi ejX + M\2eP2e.x) (j) x dxdt + JJ {Mu £ p le ^ x + Mx2eP2e,x) 4> x dxdt 
\u\>8, \v\>8 |«|<<5 

+ej I ' pi e , x <j) x dxdt + J J (M ne p le ,x + M 12e p2e,x) 4> x dxdt 

Qt |"|><5, \v\<8 

(2.37) 

and take the limit S — > extracting a proper diagonal subsequence e — > as 
follows. Using (|2.34l) one has 



lim | J J (MiiePle,* + Mi2 £ P2e, X ) 4> x dxdt ~ J J Ji<j) x dxdt\ < 
\u\>5, \v\>5 P 

lim | J J (M llE p le ,x + M 12e p2e,x) 4>xdxdt - J J Ji4> x dxdt\ + 

\u\>6, \v\>8 \v\>6, \v\>8 

+ lim | J J J\(f) x dxdt — J J Jx<f) x dxdt\ = 0, 

|«|>5, \v\>5 P 

where in the last line we used that Ji<fi x is a bounded continuous function in P. 
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Furthermore, from (12.351) one obtains 

(MuePle,x + M 12e P2e,x) <j) x dxdt 



lim 

8^0 



v=0, \u\>0 



■^H 3 pi,x + ^~M 2w I "■■■ iinl! \ •'• 



lim \ll (M lle p u . x + Mi2eP2e,x) <t> x dxdt 

5— >Q 



h. 



h. 



lim 

5^0 



\v\<5, \u\>5 



-^\ u \ 3 pi,x + ^l"! 2 ^ ) "•■ ,irdl \ ' 



v=0, \u\>0 



lim 



(5>|u|>0 



lim C 

<5^0 



I \<5>M>0 
/ 



u\ dxdt j j \u\ p l x dxdt 

5>\v\>0 



\u\ 4 dxdt J J w 2 dxdt 

\5>H>0 8>\v\>0 



< lim C 



ldxdt = 0, 



\5>|i>|>0 



where in the last inequality we used (j2.19[) and (|2.30[) . Therefore, the last two 
estimates together with f|2 .36[) imply that there exists a subsequence e — > such 
that (|2~3Tj) converge to ([2~22l . 

Similarly, defining for a fixed S > the set 72,5 = {M < 5 < one can 
estimate 



(M 2 lePle,x + M 2 2eP2e,x)4'xdxdt < CS, 



(2.38) 



\v\<8 
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and 



M 2 lePle,x + M 2 2eP2e,x ~ ^|v e | 3 p2e,ai j ",■«><// 



< C I / / ( Ml XeV {^ + 2M 21e ( M 22e - ) p lE>x p 2e 



< C 



-\u e \\v e \ 2 (Mi le p\ e x + 2M 2 i E pu. x p2e,x + M 22e p\ ex ) dxdt 



< CS 1/2 . (2.39) 

Moreover, again from (|2.30[) it follows that P2s,x P2,x weakly in L 2 (I 2 s) as 
e — » and therefore, we deduce that 

1\V E \ 3 P2S,X J <Mx£ft / / ( i|w| 3 p 2 ^ ) (2.40) 



,3 

as s — >• 0. Taking the limit 5—^0 and the corresponding diagonal sequence 
e^Oin (|2.32[) (as was done before for (|2.31[1 ) shows that it converges to (|2.23j) 
in view of (|2~38l) - (j2~40| . 

To prove (|2.26j) notice that from uo, s — > Uq, uo, e — > vq in i? 1 (SI) and (|2.7|l 
we get 



limsup / (au 2 x (x,t) + (u x + v x ) 2 (x,t)) dx < j (crul x + (u Q ^ + va^) 2 )dx. 

Q 

Since also 



u x (.,t) -> u ,a: and u-cG,*) -> «o,k weakly in L 2 (f2) 

as £ — > 0, the assertion (|2.26[) follows. 

Finally, the energy inequality (|2 . 28[) follows from (|2.7p . (|2 .29[) and the stan- 
dard weakly lower semicontinuity argument. The proof of the theorem is com- 
plete. □ 



3 Nonnegativity of solutions 

In this section we prove that the global weak solutions constructed in the pre- 
vious section are nonnegative provided the initial data uq and vq are nonnega- 
tive. Furthermore, for the system (12.11) considered with intermolecular poten- 
tials IIi(u) and II 2 (w) as in (| 1 . 3|) — (| 1 .4[) we show existence of positive smooth 
solutions. 
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3.1 Nonnegativity in the absence of intermolecular forces 

Following ideas of [3] , we define a suitable entropy in order to show nonnegativity 
of the weak solutions u and v from the Theorem 2.1, provided (|1.8p holds. 
For n £ {2, 3} we set 

A A 

9e,n( s ) = ~ j ( | r |„ + r £) i /2 ' G ^«( S ) =-j9s,n(r)dr (3.1) 

with a constant A such that A > max{|u s |, \v E \} for all sufficiently small e. 
Then one has 

G'e.nis) = 9e,n(s), G" (s) = 9 s Js) 



(|*| n + E)Va - 
Also, 

<?e,„(s)<0, G £ , n (s) >0if s< A 

and 

G e , n (s) < G (s) for all s <= K 1 (3.2) 
where Go,n = lim e _>o such that for < s < A 

f (vl — s - slog (7)), n = 2 

Go(s) = i 2(V2 + ^-Vi), » = 3 • (3 ' 3) 



Since the structure of (|1.2|) is not symmetric with respect to u and i> we use 
two entropies: one depending on u and the other on v. For a fixed 8 > one 
has 

-| / (G ei 3(v e ) + Gf e , a (w e ))da: = / (G^ i3 (u e )u c , t + G' et2 (v e )v e , t ) dx 
n n 

= / (G" 3 (We)u £ ^((Mii e +e)pi £ ^ + Mi 2£ p 2e ,x)) dx 

n 

+ / (G(,' e (u e )v eia! (M2l e pl ei a: + (M 22 e + s)p2e,x)) dx 



< S - I ((G'l 3 (u £ )) 2 (M U e + e)((M U e + e) P \ £tX + 2M 12ePlStX p 2e 



n 

+ 



(M 22e + e)p\ ex )) dx+ 5 -J ((G': e (v e )) 2 (M 22e + e)((M lle + e)p 2 u>x 

a 

+ {M 22E + e)p 2 2e . x )) dx + ^J (ul x + vl x ) dx. (3.4) 
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By definition fl3.1[) it follows that 

[G'^{u e )) 2 {M lle +e)<C, and (G'l 2 (v £ )) 2 (M 22e + e) < C, (3.5) 

where in the last inequality we used the fact that there exists a constant C such 
that M 2 2, E < C|u £ | 2 holds. Combining p.4[) with the energy inequality (|2.7I) 
and taking § < 1 one obtains 

-| | (G e , 3 (u e ) + G £)2 (t; e ))^+i^( Us , Ue ) (3.6) 
+ (1-5) J (Mus + e) P 2 U x + 2M 12ePle + (M 22e + E )j& ) (3.7) 

- ^ / + dx -^sJ ^ au ^ x + ^ Ue ' x + Ve ' x)2 ^ dx - (3 ' 8) 

n si 
This implies using Gronwall inequality 

(G e ,a( u e) + G £ ,2( v e)) d X < exp (Vj (G e ,3(uQe) + G e>2 (v 0e )) dx + -E(u 0s ,V 0e ). 



On the other hand by (jTSjl . ([23]) and (|3l2 |l -(|3T3 [l one has 

(G £i3 (M e) + G £ , 2 ("U e)) ^ < / (Go : 3("Oe) + G , 2 (v0e)) dx 



< / (G Q ,M + GoM)dx<C. 
h 

Therefore, the last two estimates imply that for all t < T 

(G e , 3 (u e ) + G e , 2 (v £ ))dx<C. (3.9) 



Finally, we prove the nonnegativity of u and v by contradiction. Assume 
there is a point (xo,io) € Qt such that u(xo,to) < 0. Since u £ —¥ u uniformly 
there exist 7 > 0, £0 > such that 

u £ (x, t ) < —7 if \x — x \ < 7, x G fi, £ < £o- 

For such x 

^ 
G e , 3 (ue(x,f )) = - J 9eA s )ds>- J g e ,3{s)ds 

u e (x,t ) —7 



u 

/ go:i( s )ds as £ — > 
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by monotone convergence theorem, where <?o,n(s) = limg^o 9e,n{s)- Since by 

JO 



Sen (s) = -oo if s < 0, n > 2 (3.10) 

it follows that 

lim G e 3 (u e (x,to)) = oo , 

which is a contradiction to (|3.9p . A completely analogous argument shows v > 
using (O and (|3TTU|) with n = 2. 

3.2 The case including intermolecular forces 

In this section we consider the system (|1.2I) - (|1.3[) in the presence of the inter- 
molecular forces given as in (|1.4[) considered with f| 1 . T[) and the initial data 
satisfying 

' [U 1 (u ) + U 2 {vo)]dx<C ll (3.11) 



<> 



where by definition 

oo 

U k (s) =- J U k (r)dT. 



Theorem 3.1. Assume that < n < m and m > 3 in Jj.^[ ). Then a positive 
smooth solution to il.2\) - (L4\ ) coupled with {1. 7\ ), iS.ll]) exists for allt G (0,T). 

Proof. Taking a suitable Holder continuous regularisation of potentials Ilfe(s) 
and proceeding as in the section 2 one can show existence of regularised solutions 
u e and v £ to (j2~2"|) considered now with (|1.2 p -(|1.4 p and (|1.7p . (|2.4p that satisfy the 
regularity properties as before. Note that the energy functional (|2.6p transforms 
in this case to 

E(u e ,v e ) = J [aul x + {u^ x + v^ x ) 2 +2U 1 {u £ ) + 2U 2 {u £ )]dx 
h 

for which the energy inequality (|2.7p still holds. Therefore, using the fact that 
Ufc(s), k = 1,2 are bounded from below, and hence 



[U 1 (u e ) + U 2 {v £ )]dx<C (3.12) 
together with (|3.11l) imply again the estimates (|2.9[) - (|2.11[) . 



£2 
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We show additionally that there exists a constant S independent of e such 
that 

u £ > 5 > 0, v s > S > hold in Q T . (3.13) 

Then proceeding to the limit e — > as in Theorem 2.1 the smoothness of the 
positive limits u and v will follow from the uniform parabolic theory and (|3 . 1 3|) 
and the statement of the theorem will be shown. Indeed, observe from (|2 . 7[) 
that 

sup / (C/^^C-,*)) + t/ 2 ( We (-,i)))rfx < C (3.14) 
te(o,T) J 

Since U2 is bounded from below one has also 

sup / Ui(u £ (-,t))dx < C. 
te(o,T) J 

Let w e (xo, t) = minn u £ (-,t). By Holder continuity of u £ we get 

u e (x,t) < u £ (x Q ,t) + C\x - xol 1 ^ 2 . 

Analogously to the proof for the single layer equation (jl.ip in [3] one obtains 
for < n < m 



C> / U 1 (u e (;,t))dx>C 2 7](u e (x ,t)) + C 3 , 



S2 



where 77(3) = — logs for m = 3, fi(s) = s 3 m for m > 3. Hence u e (a;o,t) > 
holds for all t € (0,T). The same argument works for minn v £ (-, t). Therefore 
(j3~TI)) is true. □ 

4 Existence of nonnegative weak solutions in the 
Navier-slip case 

In this section we show that solutions u £ and v £ to the regularised system (|2.2p - 
(12. 4[) considered now with the Navier-slip mobility matrix 



M 



\u\ 2 + e \u\\v\ 
\u\\v\ (l + a)|7;| 2 +e; 



converge to global nonnegative weak solutions to (|1.2I) considered with (|1.6[) , 
(|2.ip and (|1.7|) - (|1.8|) . Note that the case when intermolecular forces are present 
i.e. for (jl.2p considered with (jl.3l) - (jl.4l) and (jl.6p proceeds then exactly as in 
Theorem 3.2 for the no-slip case. 



1G 



The dissipation (|2.7|l of the energy functional (|2.6p and the corresponding a 
priori estimates (|2.9|) — (|2.13|) and f|2 . 16[) stay true in the Navier-slip case as well. 
Therefore, (|2.18l) holds again up to a subsequence as e — > 0. The following the- 
orems that thus obtained limits u and v are nonnegative global weak solutions. 



Theorem 4.1. Functions (u,v) satisfy for any T > the following properties: 

u, V eC^ 1/8 (Q T ), u,veC 4 x :l(P), (4.1) 
M u p M + M 12 p 2 ,x, M 21 pi, a + MaaPa,* G i 2 (P), (4.2) 
MV* e L 2 (i?), |v| 2 |>2,* G L 2 (5); (4.3) 



Qt 



where_P = Qt\({u = 0} U {v = 0} U {t = 0}), i? = QtC\{v — 0}n{|u| > 0} 
S = Qt H {u = 0} n > 0}. Furthermore, there exist functions w% G L 2 (R) 
and W2 G L 2 (S), such that (u, v) satisfies U.fy) in the following sense: 

J J ' v4t+ J J (Mnpi, x + Afi 2 p2 !2; ) 0x (4.4) 
(|«| 2 Pi,x + |w|wi) 0s = 0, 

■'"/ + / / i -U 2 ipi,x + M 22 p 2 ,x) 4>x+ J J (\v\w 2 + (1 + a)\v\ 2 p 2 ^)<j) x = 

s 

(4.5) 



Qt 



for all 4> G Lip(Qt), 4» = near 4 = and i = T; 

u(x, 0) = uq(x), v(x, 0) = Vo(x), x G ft, (4-6) 

IK',*)IUi(n) = H u o||li(o), \\v(-,t)\\ L i m = \\u \\ L i m , (4.7) 

u x (-,t) —> uq.x and v x (-,t) —> Vq iX strongly in L 2 (ft) as t — > 0, (4-8) 



and 



u and v satisfy \1. 7| ) at all points of the lateral 
boundary, where and v 0. 



(4.9) 



Finally, the following energy inequality holds 

E(u(;T),v(;T))+ ff (M 2lP l x + M 22 p\ x ) <E(u ,v ). (4.10) 



Proof. The assertions (|4.1[) - (|4.2[) and (|4.6[) -(|4.8 j) follow exactly as in the proof 
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of Theorem 2.1. Using (|2.10p one observes 

J I {M lle p\ ex + M 22e p\ ex )dxdt < C - 2 J J Ml2ePls, X P2 e ,xdxdt 
Qt Qt 

= C-2jj{\u e \\v e \ Px ^)d X dt 

Qt 

By Young's inequality 

2\u e \\v e \p lEtX p 2s ,x) < -^-^\Ue\ 2 ple, x + ( X + f ) \ V ^vle,x 

and hence one obtains 

J \u s \ 2 p\e,xdx < C, J \v e \ 2 pl,xdx < C, (4.11) 

Qt Qt 

J j {M lle p\^ x + M 22e p 2 2e , x )dxdt < C. 

Qt 

From the last inequality (|4.10j) then follows. 

Next, for <f> as in (|4T5|) -(|43 |) one writes again f|2.31|) — ()2.32[) . Considering the 
set li s and I 2} s as in the proof of Theorem 2.1 and using the estimates (|4. 1 1[) 
one can show in the analogous manner that (12.31[) and (|2.32l) converge up to a 
diagonal subsequence as S — > and e — > to (|4.4p and (|4.5p respectively. □ 

Theorem 4.2. The global weak solutions u and v constructed in the previous 
theorem are nonnegative. 

Proof. The proof proceeds similarly to the argument for nonnegativity of the 
weak solutions in the no-slip case presented in the section 3.1. The following 
estimates 

(G'l 2 (u £ )) 2 (M ll£ +e)<C, and (G^ 2 (« £ )) 2 {M 22e + e) < C, 

and 

d f 1 d 

(G e , a (u e ) + G e , 2 (v e ))dx + -— E{u £ ,v e ) 



dt J v ' v ' ' v " 2dt 



+ {l-5) J (M ll£ +e)p\ E , x + 2M l2s p l£ ,xP2s,x + (M 22e + e)p\ e x ) dx 

-hi + v ^ dx ' 
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are used as analogs to (|3.5[l - (|3.6[l in this case in order to obtain the crucial 
estimate 

[ {G e , 2 (u e ) + G E!2 (v £ ))dx< C. 
h 

The rest of the proof proceeds proceeds exactly as in the last paragraph of the 
section 3.1 but now using (|3.10[) with n = 2 for both functions u and v. □ 

5 Conclusion and discussions 

In this article we showed existence of nonnegative global weak solutions for the 
coupled lubrication systems corresponding to the cases of no-slip and Navier-slip 
conditions at both liquid-liquid and liquid-solid interfaces. Our results can be 
generalised in a straight forward way to the system (ll.2[) considered with the 
mobility matrix 

^ / |u 3 + biu 2 \ u2y + b\uv \ 

M = - . (5.1) 

M \\u 2 v + b x uv f v 3 +uv 2 + b x v 2 + b(fi + l)v 2 J 

corresponding to the weak-slip conditions at the both interfaces. As it was 
shown recently in [T3] the latter model incorporates both the no-slip and the 
Navier-slip models (|1.5[) . (|1.6[l as limiting cases as the slip lengths 6, b\ tend 
simultaneously to zero or infinity, respectively. 

One needs to point out that we obtained a slight difference between the 
weak formulations in the no-slip and Navier-slip cases (compare Theorems 2.1 
and 4.1). Due to the fact that Mn and M 2 2 components in (11.61) depend only 
on u or v, respectively, in contrast to the no-slip case (|1.5[) there is no an analog 
of estimate (|2.38[) in the Navier-slip case. Therefore, an additional (so far not 
identified in terms of solutions u and v) function u>2 appears on the singular 
set S in the latter case. The same problem persists also in the weak-slip case 
because the leading orders of the mobility matrix components on the set S 
coincide with those for the Navier-slip case. 

In this sense we have "more regularity" for the weak solutions in the no-slip 
case then for ones the in weak- or Navier-slip cases. This interesting observation 
should be understood better in future in view of the fact that for the single 
lubrication equation (jl.ljl the no-slip case is known to be more singular from 
both physical and analytical points of view then the slip cases. At the same time 
we've become aware of an alternative proof for the existence of weak solutions in 
the no-slip case in |12j for which the authors have shown the same regularity as 
we in the Navier-slip case. But also in our weak formulation for the no-slip case 
remains an open question weather not yet identified in terms of the solutions 
function w\ vanishes on the singular set R. Another observation appearing as 
well due to different component structures of the mobility matrices (( 1 . 5|) and 
(|1.6[) is that we have stronger entropy for u then for v in the case of the former 
matrix whereas the entropies are the same for the latter one. 



19 



Additionally in contrast to the existing results for the single lubrication 
equation (jl.ll) (see e.g. [3 [5]) we are not aware if the constructed weak solutions 
for the systems (ll.2p should necessarily posses zero contact angles. This is due to 
an absence so far of the strict entropy dissipation inequality for the two layered 
systems p.2[) which was shown before to hold for (jl.ip . Combined energy- 
entropy dissipation inequalities derived here (as e.g. (I3.6[l ) do not imply H 2 a 
priori estimates on the solutions. 
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